Anomalous Scaling in the Anisotropic Sectors 
of the Kriachnan Model of Passive Scalar Advection 
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Kraichnan's model of passive scalar advection in which the driving (Gaussian) velocity field has 
fast temporal decorrelation is studied as a case model for understanding the anomalous scaling be- 
havior in the anisotropic sectors of turbulent fields. We show here that the solutions of the Kraichnan 
equation for the n order correlation functions foliate into sectors that are classified by the irreducible 
representations of the SO(d) symmetry group. We find a discrete spectrum of universal anomalous 
exponents, with a different exponent characterizing the scaling behavior in every sector. Generi- 
cally the correlation functions and structure functions appear as sums over all these contributions, 
with non-universal amplitudes which are determined by the anisotropic boundary conditions. The 
isotropic sector is always characterized by the smallest exponent, and therefore for sufficiently small 
scales local isotropy is always restored. The calculation of the anomalous exponents is done in two 
complementary ways. In the first they are obtained from the analysis of the correlation functions of 
gradient fields. The theory of these functions involves the control of logarithmic divergences which 
translate into anomalous scaling with the ratio of the inner and the outer scales appearing in the 
final result. In the second way we compute the exponents from the zero modes of the Kraichnan 
equation for the correlation functions of the scalar field itself. In this case the renormalization scale 
is the outer scale. The two approaches lead to the same scaling exponents for the same statistical 
objects, illuminating the relative role of the outer and inner scales as renormalization scales. In 
addition we derive exact fusion rules which govern the small scale asymptotics of the correlation 
functions in all the sectors of the symmetry group and in all dimensions. 



I. INTRODUCTION 

The aim of this paper is twofold. On the one hand we 
are interested in the effects of anisotropy on the universal 
aspects of scaling behavior in turbulent systems. To this 
aim we present below a theory of the anomalous scal- 
ing of the Kraichnan model of turbulent advection [Q in 
anisotropic sectors which are classified by the irreducible 
representations of the SO(d) symmetry group. On the 
other hand we are interested in clarifying the relationship 
between ultraviolet and infrared anomalies in turbulent 
systems. Again, it turns out that the Kraichnan model is 
an excellent case model in which this relationship can be 
exposed with complete clarity. As is well known by now, 
the Kraichnan model describes the advection of a pas- 
sive scalar by a velocity field that is random, Gaussian 
and delta-correlated in time. The correlation functions of 
the field scale in their spatial dependence, and the main 
question is what are the scaling (or homogeneity) expo- 
nents of the statistical objects of the scalar field that are 
induced by the given value of the scaling exponent of the 
advecting velocity field. 

The two issues discussed in this paper have an im- 
portance that transcends the particular example that 
we treat in detail in this paper. The first is the role 
of anisotropy in the observed scaling properties in tur- 
bulence. We have shown recently that in the presence 
of anisotropic effects (which are ubiquitous in realistic 
turbulent systems) one needs to carefully disentangle 



the various universal scaling contributions. Even at the 
largest available Reynolds numbers the observed scaling 
behavior is not simple, being composed of several contri- 
butions with different scaling exponents. The statistical 
objects like structure functions and correlations functions 
are characterized by one scaling (or homogeneity) expo- 
nent only in the idealized case of full isotropy, or infinite 
Reynolds numbers when the scaling regime is of infinite 
extent. Anisotropy results in mixing in various contri- 
butions to the statistical objects, each of which is char- 
acterized by one universal exponent, but the total is a 
sum of such contributions which appears not to "scale" 
in standard log-log plots. By realizing that the correla- 
tion functions have natural projections on the irreducible 
representations of the SO (3) symmetry group we could 
offer methods of data analysis that allow one to measure 
the universal scaling exponents in each sector separately. 
In this paper we show that this foliation is an exact prop- 
erty of the statistical objects that arise in the context of 
the Kraichnan model. 

The second issue that transcends the particular exam- 
ple of the Kraichnan model is the identification of the 
renormalization scales that are associated with anoma- 
lous exponents. As has been already explained before, 
the renormalization scale which appears in the correla- 
tion functions of the passively advected scalar field is the 
outer scale of turbulence L. On the other hand, the the- 
ory of correlations of gradients of the field expose the 
inner (or dissipative) scale a an additional renormaliza- 
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tion scale. Having below a theory of anomalous scaling 
in various sectors of the symmetry group allows us to 
explain clearly the relationship between the two renor- 
malization scales and the anomalous exponents that are 
implied by their existence. Since we expect that Kol- 
mogorov type theories, which assume that no renormal- 
ization scale appears in the theory, are generally invali- 
dated by the appearance of both the outer and the inner 
scales as renormalization scales, the clarification of the 
relation between the two is important also for other cases 
of turbulent statistics. 

The central quantitative result of this paper is the ex- 
pression for the scaling exponent £n which is associated 
with the scaling behavior of the n-order correlation func- 
tion (or structure function) of the scalar field in the Pth 
sector of the symmetry group. In other words, this is the 
scaling exponent of the projection of the correlation func- 
tion on the £'th irreducible representation of the SO(d) 
symmetry group, with n and £ taking on even values only, 
?i = 0,2,... and £ = 0,2,...; 



C 



n(n + d) (d + !)£(£ + d- 2) 



2(d + 2) 2(d + 2)(rf-l) 



+ 0(e 2 ). 
(1.1) 



The result is valid for any even £ < n, and to O(e) where 
e is the scaling exponent of the eddy diffusivity in the 
Kraichnan model (and see below for details). In the 
isotropic sector {£ = 0) we recover the well known re- 
sult of It is noteworthy that for higher values of £ 
the discrete spectrum is a strictly increasing function of 
£. This is important, since it shows that for diminishing 
scales the higher order scaling exponents become irrele- 
vant, and for sufficiently small scales only the isotropic 
contribution survives. As the scaling exponent appear in 
power laws of the type (i?/A)'=, with L being some typ- 
ical outer scale and R <C A, the larger is the exponent, 
the faster is the decay of the contribution as the scale 
R diminishes. This is precisely how the isotropization 
of the small scales takes place, and the higher order ex- 
ponents describe the rate of isotropization. Nevertheless 
for intermediate scales or for finite values of the Reynolds 
and Peclet numbers the lower lying scaling exponents will 
appear in measured quantities, and understanding their 
role and disentangling the various contributions cannot 
be avoided. 

The organization of this paper is as follows: In Sect. 2 
we recall the Kraichnan model of passive scalar advec- 
tion, and introduce the statistical objects of interests. In 
Sect. 3 we set up the calculation of the correlation func- 
tions of gradients of the field. It turns out that it is 
most straightforward to compute the fully fused correla- 
tion functions of gradient field, as these objects depend 
only on the ratio of the outer and inner scales. We com- 
pute these quantities and their exponents to first order 
in e. We introduce the appropriate irreducible represen- 
tations of the SO(d) symmetry group and evaluate the 
scaling exponents in all its sectors. In Sect. 4 we turn to 



the correlations functions of the passive scalar field itself, 
and compute the scaling exponents of the structure func- 
tions in the presence of anisotropy, again correct to first 
order in e. To this aim we compute the zero modes in all 
the sectors of the symmetry group. One of the interesting 
points of this paper is that the results of this calculation 
and the calculation via the correlations of the gradient 
fields gives the same results for the scaling exponents if 
one accepts the fusion rules. To clarify the issue we prove 
the fusion rules here in all the sectors of the symmetry 
group by a direct computation of the fusion of the zero 
modes. In Sect. 5 we offer a summary and a discussion. 



II. KRAICHNAN'S MODEL OF TURBULENT 
ADVECTION AND THE STATISTICAL OBJECTS 

The model of passive scalar advection with rapidly 
decorrelating velocity field was introduced by R.H. 
Kraichnan Q already in 1968. In recent years it 
was shown to be a fruitful case model for understanding 
multi-scaling in the statistical description of turbulent 
fields. The basic dynamical equation in this model is for 
a scalar field T(r, t) advected by a random velocity field 
u(r,t): 



[d t - Ko V 2 + u(r, t) ■ V] T(r, t) = f(r, t) 



(2.1) 



In this equation f{r,t) is the forcing. In Kraichnan's 
model the advecting field u(r,t) as well as the forcing 
field f(r,t) are taken to be Gaussian, time and space 
homogeneous, and delta-correlated in time: 



f(r,t)f(r>,t>) = $(r-r')6(t-t'), 
(u a (r, t)u p {r', t')) = W af3 (r - r')5{t - t') , 



(2.2) 
(2.3) 



Here the symbols 7TT and (■ • •} stand for independent 
ensemble averages with respect to the statistics of / and 
u which are given a priori. We will study this model 
in the limit of large Peclet (Pe) number, Pe= UjyA/K , 
where Ua is the typical size of the velocity fluctuations on 
the outer scale A of the velocity field. We stress that the 
forcing is not assumed isotropic, and actually the main 
goal of this paper is to study the statistic of the scalar 
field under anisotropic forcing. 

The correlation function of the advecting velocity 
needs further discussion. It is customary to introduce 
W af3 (R) via its fc-representation: 



eD 



d d p 



O d J p- 
A- 1 



+ c 



p a ^{p) 



Sal: 



P a/3 (p) exp(-ip- R) , (2.4) 



(2.5) 



were P a ^(p) is the transversal projector, fid — (d — 
l)il(d)/d and £l(d) is the volume of the sphere in d di - 
mensions (i.e. 0(2) = 2n, 0(3) = 47r). Equation (2.4) 
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introduces the four important parameters that determine 
the statistics of the driving velocity field: A and A are 
the outer and inner scales of the driving velocity field 
respectively. The scaling exponent e characterizes the 
correlation functions of the velocity field, lying in the in- 
terval (0, 2). The factor D is related to the correlation 



(2.6) 



function (2.3) as follows: 

W*f , (0)=D5 af} (A e 

The most important property of the driving velocity field 
from the point of view of the scaling properties of the 
passive scalar is the tensor of "eddy diffusivity" Q 

K f(R) = 2[W Q/3 (0) - W aP (R)} . (2.7) 

The scaling properties of the scalar depend sensitively 
on the scaling exponent e that characterizes the R de- 
pendence of k"P(R): 



K f(R) ex [A e - A £ ]^ 



(R) cx R* 



for R > A , 

R a R !3 i 







R 2 



(2.8) 
A < R < A . 



-1 + e 

We are interested in the scaling properties of the scalar 
field. By this we mean the power laws characterizing 
the R dependence of the various correlation and response 
functions of T(r,t) and its gradients. We will focus on 
three types of quantities: 

1) "Unfused" structure functions arc defined as 

S n ( ri ,r u ...r n ,r n ) = ([T(r 1; t) - T(r x , t)} (2.9) 

x [T(r 2 ,t) - T(r 2 ,t)} . . . [T(r n ,t) - T(r n , <)]) , 

and in particular the standard structure functions are 

S n (R) = ([T(r + Rt)- T(r, t)] n ) . (2.10) 

In writing this definition we used the fact that the sta- 
tionary and space-homogeneous statistics of the velocity 
and the forcing fields lead to a stationary and space ho- 
mogeneous ensemble of the scalar T. If the statistics is 
also isotropic, then S n becomes a function of R only, in- 
dependent of the direction of R. The "isotropic scaling 
exponents" Q n of the structure functions 

S n (R)(xR^, (2.11) 

characterize their R dependence in the limit of large Pe, 
when R is in the "inertial" interval of scales. This range 
is A, n <C R <C A, L where n is the dissipative scale of the 
scalar field, 



(2.12) 



2) In addition to structure functions we are also inter- 
ested in the simultaneous nth order correlation functions 
of the temperature field which is time independent in 
stationary statistics: 

^ n ({r m }) = (T(r 1 ,t)T(r 2 ,t)...T(r n ,t)) , (2.13) 

were we used the shorthand notation {r m } for the whole 
set of arguments of nth order correlation function Tm 
n,r 2 . . .r„. 



3) Finally, we are interested in correlation functions of 
the gradient field VT. There can be a number of these, 
and we denote 



n 

H^{{r m }) = (j{^T{r u t) 



i=l 



(2.14) 



where {a m } is a set of even n vector indices {a m } = 
ai, a.% . . . a. n . We introduce also one-point correlations 
which in the space homogeneous case is independent of 
the space coordinates: 



H^^HZ ia2 - a »({r m = r}). 



(2.15) 



The tensor U^ m} {{r m }) can be contracted in various 
ways. For example, binary contractions oi\ = a 2 ,ct3 = 
«4, etc. with 7*i = r 2 , = etc. produces the corre- 
lation functions of dissipation field |VT| 2 . 

When the ensemble is not isotropic we need to take into 
account the angular dependence of R, and the scaling 
behavior consists of multiple contributions arising from 
anisotropic effects. The formalism to describe this is set 
up in Appendix A and in the forthcoming Sections. 

The correlation functions T n satisfy equation Jl[ 



= 1 E ^ri-^fKWy), ( 2 - 16 ) 

{i#j}=l 

where {r m } m ^ij is the set off all r m with m from 1 to 
n, ex cept of m = i and m = j. Substituting k"P(t) from 
Eqs. (§J, |3) one gets: 



-«E V ?+ E W Q V*-^)V?vf].F n ({r m }) 

n 

= 3 E ^-r 3 )T n - 2 ({r m } m ^), (2.17) 



where 



(2.18) 



Here we used that in space homogeneous case X)"=i = 
and therefore 

|E V *| =E V '+ E v?vf = o. 

«=i i=i {i^-}=i 

Consider the fc-Fourier transform of T n which is defined 
as: 



(2n) d s(j2k s )F n ({k m }) 



s=l 



(2.19) 



[ ]J rfr mexp(ifc m -r m ) f„({r ra }) 

m— 1 



3 



Here the 5-function app lies t o a homogeneous ensemble in which J r n ({r m }) depends only on differences of coordi- 
nates. For F n ({k m }) Eq. fl2.17|) yields: 



a- 1 



K F n ({k m })^2k^ + ^ f ^P aP (p) £ k?kfF n (k i +p,k j -p ) {k m } m7ti , j ) = $n({k m }) 
d A -i P {*7*3}=1 



d n 



$n({fem}) = X] + fej)i ? n-2({fem}m#i,j) , for 71 > 2 

{i#j}=l 

#(fe) = / dRexp(ik ■ R)<f>(R) 



Here l>(fc) is the Fourier transform of $(-??) and $2(fe) = $(fe). Equation ( 2.20| ) may be rewritten as: 



a- 1 



A- 1 



Fn(ki + P, fcj P, {fcm}m^i,j) ' 



This equation will serve as the basis for future analysis in Sect. III. 



(2.20) 

(2.21) 
(2.22) 

(2.23) 



III. SCALING OF THE TEMPERATURE GRADIENT FIELDS 



A. Basic Equations in k-representation 



It appears that Eq. ( |2.23 ) is as dif ficult to solve as Eq. ( 2.1 6| ) . In fact, very important information about scaling 
behavior may be extracted from Eq. (2.23) for small e 8. In order to develop our method we will analyze first the 
simultaneous, n-point correlation functions of the gradient fields 'H-if"'^ ({r m }) and Hn"" 1 ^ of Eqs.(2.14,2.15): These 
objects are expressed in terms of F n ({k m }) as follows: 



/n n 
II [d d kik ai exp(tfc i .r i )]i ? „({fcm})«(53 fc «)' 
i=l s=l 
„ n n 



(3.1) 
(3.2) 



From this and Eq. (2.23) one gets 



rU n h a 'd d h / ™ \ r ri d n P al3 (v) Y!-^ i , k a k 

s- 1 A-i S 



^ n ' x n 



K=iK°d d k s $„({fc m })^^ 



(I» 



(2 1 r)(»-'>" «e;,.*2 ._, 

Shifting the dummy variables fe$ — p — ► ki and fej + p — > fcj we have another representation of this equation 

d d p (fcf -p Q -)(^ + P ^)P afi {p)kfk l3 ] 



}m?ij) (3-3) 
(3.4) 



/sJW (27r)C»-i)d V ; .^_J 2p 2 + 2p-(fc,-fc l ) + ELifc s 2 



s—l,s^i,j 



(3.5) 



In order to analyze this equation further we choose to nondimensionalize all the wave vectors by A. We write k s = Ak s , 
p = Ap etc, and for simplicity drop the tilde signs at the end. We simplify the appearance of the equation further by 
introducing the definition of the dimensionless function 
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Ap t fi i M m teij*i>*J>P)- J ^ 2 p2 + 2p.(fc,-fc l )+E:=ifc s 2 ' 
where dp stands for integrating over all the angles of the unite vector p = p/p. The resulting equation is 

A/A 

J 1 j »=1 1 ^ 

n 

x II K 3 Fn({k m }) + *[ a ^ , (3.7) 



where the dimensionless factor g is 
3 = k + D(A< - A £ ) ' 



(3. 



In fact, we should recognize that the natural expansion 
parameter is actually not g, but g, where 



9 



A/A 



P 



edp 
i+i 



Evaluating the integral we find 
L>[A e - A e ] 
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(3.9) 



(3.10) 



K + £>(A e - A e ) 

We will see below that g can take on very different values 
in different limiting cases. In particular it can be of 0(e) 
or of 0(1) depending on the order of limits. The relevant 
limit for the physics at hand will be discussed below. At 
this point we perform a calculation of Hn° m ^ to first or- 
der in g in all the sectors of the symmetry group. 



B. Theory to first order in 



F n ({r m }) = £ F n , q ({r m }) , = £ . 

g=0 q=0 

(3.11) 

Here F n „ is the result of the g-step of the iteration proce- 



dure of Eq. (2.23), F Jhq oc g q . There are two contributions 



for each term of order g q . One arises fro m su bstituting 
F n ,q-i into the integral on the RHS of ( 2.23| ), and the 
sec ond a rises from F n _2,q which appears in according 
to (2.21). Correspondingly, also H n a has two contribu- 
tions. One is obtained from Eq. (3.7) when we substitute 
F n .q-i in the integral on the RHS, and the second when 
we substitute F n _2,q in the term denoted as \& n . 

In this Sect, we compute explicitly H n> i. Analyzing 
the relative importance of these two contributions to H n \ 
we found that the 2nd contribution is negligible compared 
to the 1st when A -C L. In other words, we can disregard 
the contribution to H Ui i which arises from F n _2,q- This 
means that for the sake of our iteration procedure we 
can replace $„ in ( 2.23j ) by the quantity <j> n ,o which is of 



0(g ). This means that instead of Eq. (2.23) we iterate 



The theory for F n and H n can be formulated itera- 
tively, resulting in the following series: 



D A / d d T n ^k 13 

F n ({k m \) = -i±l J _P p^(p) ^ {l ^)r\l 3 F n (ki + p, fc, - p, {fc m W»,j) + Fn,o({fcm» , (3.12) 
A" 1 1 

F n , ({k m }) = $ "^ fcm >) , $„ >0 ({fe m }) = ^1 V $(fc*)^ + fei)^-2,o({fe m }m^,i) • (3.13) 



Thus we are interested in calculating 



A/A 

1 j s=l {«#}=! i ^ s=l,s#i,j 



.}) • (3-14) 



Recall that the function $(fe) is constant for kL <C 1 and it vanishes for fcL 3> 1. Therefore the leading contribution 
to the integrals (3.14) over ki come s from the region fcjL < 1. In integral (3.14) p > 1/A and in all our approach we 
consider L ^> A. Therefore in (3.14) p S> % and this equation may be simplified up to: 
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H 



{c m } 



IE 



d d k. 



A/A 



(2tt)(" 



-l)d 



s=l 



P 



n *z 



'i^,0({fem}) 



(3.15) 



where now 



(3.16) 



is the con stant tensor that obtains from the tensor func- 
tion ( |3.6| ) when all k s <C p. Performing the all the wave- 
vector integrals we observe that the explicit e is canceled 
by integral over p. Accordingly 



Defining now A$ via the relation 



we conclude that 

w _ 2n(d + n) _ (d + 1)1(1 + d- 2) 
n d + 2 2(d-l)(d + 2) 



(3.20) 



(3.21) 



{i#j'}=l 



(3.17) 



D. Interpretation of the result 



An actual integration in ( 3.16| ) yields 

_ ^Qia^ftft (d + 1) — S ai/3i S a:j p j — 5 ai j3 j 5 aj j3 i 



A 



2(d + 2)(d- 1) 



(3.18) 



C. Analysis in all the anisotropic sectors 



We note that Eqs. ( [2.17 ) contain only isotropic oper- 
ators. On the other hand $(r^ — rj) can be anisotropic, 
depending on the direction of the vector — r j . Since 
he equations are linear, we can expand all the objects 
in terms of the irreducible representations of the SO(d) 
group of all rotations, and be guaranteed that the solu- 
tions foliate in the sense that different irreducible repre- 
sentations cannot be mixed. This considerations are valid 



for all the equations in this theory, including Eq. (3.17). 
To know which irreducible representations we need to 
use in every case one has to consult Appendix A. Af- 
ter doing so one notes that for any order q, the tensors 
Hn"^ are constant tensors, fully symmetric in all their 
indices. Using the exposition of Appendix A we know 
that the projections on the irreducible representations of 
the SO(d) symmetry group must be of the form 



rr{<x m } _ X W B {a m } 



(3.19) 



Our first order calculation is aimed at finding the ratio 



A 



Substituting ( |3 . 1 8| , |3 . 1 9| ) in ( |3.17| ) we find 



H, 



{« m } 



n ' 1 ' e 2(d + 2)(d-l 



(d+i^^s^A^ 
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(d + 2)(d-l) 



(d + l)z n j 



n(n — 1) 



A a i,a,n 



To interpret the result (|3.20H3.2l| ) we should observe 
that the nature of the theory that we develop depends 
on the order of the limits that we take. We should recog- 
nize that the quantity Hn"" 1 ^ does not exist without an 
inner (ultra-violet) cutoff. We are thus interested in lim- 
iting values of g subject to the condition that r\ is finite. 
Thus one order of limits that makes sense is A — > first 
( corresponding to the Reynolds number going to infin- 
ity first), and then e going to zero second, b ut kee ping r\ 
fixed [for example by controlling kq in Eq. (fU2j )1. An- 
other order of limits is e — > first, (still keeping rj fixed, 
but very small) with A being fixed and larger than r\. 

If we have A — > first, and then when e — > second 
we find that the expansion parameter is close to unity: 



1 



(1)' 



for e In 



(1) 



< 1 



(3.22) 



Thus we cannot stop at ( |3.20| ), and we are forced to con- 
sider higher order terms in the expansion in g and appro- 
priate resummations. This is done in Subsect. F. On the 
other hand, if e — > first, we find an apparently "small" 
expansion parameter that is proportional to e: 



■'"'I 



for 



-1 



> 1 . 



(3.23) 



E. Exponentiating using Renormalization Group 
Equations 



Using Eq. ( |3~23| ) in Eq. ( jTgCj ) we get: 



H 



n,0 ■ 



(3.24) 



If we expect that H*p is a scale invariant function of 
A/ A we can interpret Eq. (3.24) as the beginning of an 
expansion that can be re-summed into a power law 



G 



5 A« 



H 



n,0 



(3.25) 



Of course, this is hardly justified just by examining the 
0(e) term, since one can have more than one branch 
of scaling exponents proportional to e. If we have m 
branches only the analysis up to 0(e m ) can reveal this. 
We return to this issue in the next Subsection. An addi- 
tional issue is the magnitude of A that can be arbitrarily 
small, making any re-exponentiation even more dubious. 
To overcome this problem one usually invokes the renor- 
malization group equations to justify the exponentiation. 
We shortly present this method next. In doing so we 
want to argue that for the case in question there is noth- 
ing more in this approach than direct re-exponentiation 
as long as higher order in e are not included. 

Within the renormalization group method || one con- 



siders Eq. (EyJ) as the "bare" value of H$, H 



One then seeks a multiplicative renormalization group 
by defining a renormalized function 



= Z H ( Pl A,X)H^ B (A,X) . (3.26) 



Here \i is a "running length" , and the only [i dependence 
of the RHS is through the Zh function. Defining Zh so 
that it eliminates the dependence of the LHS on A and 
setting the initial condition H nt n(A, A, . . .) = H 
get: 



n,0 



Z H (fi,X) =l + eA^ ln(A//i) 
From equation ( |3.26 ) we get: 

dHH ( X (M,A,...)) 



din ji 



1h 



(3.27) 



(3.28) 



where 



1 H = -^+0(6 2 ). (3.29) 
Solving the differential equation ( |3.28 ) we get: 



<W,A,. 



(3.30) 



Exponentiating ( 3.27 ) and solving Eq. ( 3.26| ) in favor of 



H^ B we recover Eq. (3.25). Note that the inner scale in 
this case is A, since A > r/. This fact casts an additional 
doubt on this limit of the theory, since it misses alto- 
gether the existence of the Batchelor regime Q between 
A and r\. For all these reason we tend to disqualify ( 3.25| ) 
despite the relative simplicity of its derivation. We turn 
next to the other limiting procedure. 



F. Theory for A -» first 



Considering the limit A — > first, Eq. (3.20) is still 
valid, but now g is of order unity, and we cannot justify 
re-exponentiation by any stretch of the imagination. 



FIG. 1. The graphic representation of Eq. (3.3)(Psi4). 



a, a, 
k , k 



, IXi (XX. 
k i+Pj 




FIG. 2. The graphic representation of Eq. (3.8) (HA). 



The small parameter e seems to have disappeared. This 
forces us to consider all the higher order terms in g to 
understand how to resum them. We will see that at the 
end e reappears. 

The re-summation of the g dependence is aided signif- 
icantly by some graphic representations of the relevant 
equations and their perturbative solutions. In Fig. Fy 

rfaw 



we represent graphically the definition (3.2) of Hn 
in terms of F n ({k m }). This helps us to introduce the 
basic diagrammatic notations. 

A solid long rectangle with n a rn wavy lines stands 
for Hn"" 1 ' , while the elongated solid ellipse represents 
F n . The k s wave- vectors are denoted by arrows, and the 
dots represent multiplications. The vertical line connect- 
ing all the arrow heads stands for the integration with a 
delta function over the sum of all the k vectors. Fig. || is 



a graphic notation of Eq. (3.7). The little ellipse with a 
vertical arrow designated by p stands for the tensor func- 
tion A^p J ({fcsjs^i .jki, kj, p). This ellipse is involved in 
the integration over the vector p in the loop to the left of 
the ellipse with a weight consisting of the sum of squares 
of the k vectors. 

The thin elongated ellipse in the second term on 
the RHS stands for the zero'th order term F ny o, cf. 
Eq. (3.13). The actual values of the wave- vectors are 
indicated in this diagram. In later diagrams, Fig. 3, we 
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drop this obvious notation. In Fig. |^ we display the per- 
turbative solu tion which results from the iteration pro- 
cedure in Eq. ( 3.12 ), the result of which is substituted in 
Eq. Q. 




FIG. 3. The graphic representation of the iteration scheme 

(£) 

The first diagram on the RHS is H n , whereas the 



second is the first order term H^\, Eq. ( 3.20| ). The last 
diagram shown is i?„ 2 



.4' 



One should note that when the 
p.p. ellipse appears once we have a sum over all pairs 
i,j indices with i ^ j. When it appears twice there is a 
double sum, with respect to the pairs i ^ j and £ ^ m. In 
analyzing such diagrams one needs to identify three dis- 
tinct possibilities. These are denoted as case (a), i = £, 
j = m, case (b) i = £ but j ^ to, and case (c) where all 
the indices are different. There are two integrals over p 1 
and p 2 in the loops to the left of the corresponding ellipse. 
We refer to the functions A of Eq. (^^) which appears in 
these integrals as A\ and Ai respectively. Analyzing the 
integrals it is useful to separate the discussion to region 
I in which p\ > P2, and region II in which p2 > Pi- In 
region I we can neglect the contribution of p 2 and all k s 
with respect to p 1 . Accordingly we have two independent 
integrals and sums, and the result is therefore 



H 



\[gA^fH%, (region I) 



(3.31) 



where the factor 1/2 stems from the fact that the volume 
of region I is a half of the whole volume of (p 1 ,p 2 )-space. 
In region II we should distinguish between the cases (a), 
(b) and (c), for which the evaluation of A^ will be differ- 
ent. In case (c) Eq. ( ]3.6| ) shows that A 2 is of the order 
of p\jp\i which is small. In case (b) A 2 is of the order 
of P2/P1, which is still small. Only case (a), in which the 
the loops appear as two rungs on the same ladder, we 



have A2 of the order of unity. The actual calculation of 
this integral is presented in Appendix B, with the final 
result 

H% = W^H% , [region II, case (a)] . (3.32) 
Together the second order result for H n 2 is 



"n,2 — 2^ 



A w + ( A w 

71 1 \ n 



H 



n,0 



(3.33) 



Our aim is to find the fully resummed form, correct to 
all order in g and A„ , of H$ . We can express it in the 
form 



H 



m = 



K{g,A^)H% 



(3.34) 



where the function K(g, A„ ) is represented as the dou- 
ble infinite sum 

00 00 

K(g, A) = 1 + A m ]T An, s .g s ■ (3.35) 

m— 1 s—m 

Up to now we have information about -Di,i = 1 and 

D1.2 = D2.2 = \- 

In Appendix C we derive the following recurrent rela- 
tion for the higher order terms 



_ _ 1 1 ^ 

e Q £ ^ 



(3.36) 



Using ( |3~36| ) in Eq. ( |3~35|) we find the contribution pro- 
portional to A: 



00 ~ « 
9 



K l {g,A)=AY^ y - = -Ahx{l-~g) 



(3.37) 



Considering all the terms quadratic in A, and using the 
recurrent relations to determine D 2 s we find 



OO ~ „ S—l 



K2{g,A) = A^Y.\ 



s=2 3 Pi q 



(3.38) 



This double sum is computed in Appendix C with the 
result 



K2{g,A)= 1 -[-A\n{l-~g)¥ 



(3.39) 



The general result can be derived using similar techniques 
with the result 



1 



K m (g,A) = —[-A\n(l-g)Y 



(3.40) 



Accordingly we conclude with the series for K (g, A): 
-Aln(l-g)] n 



K(~g,A)=J2 



m=0 



exp[-Aln(l - g)} . 

(3.41) 
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Using now Eq. (3.22) one finds 



lu,. .[) = [- 

.V 



Finally, using Eq. ( 3.34 ) we have the final result 



H 



(i) 

ra.O 



(3.42) 



(3.43) 



We are pleased to find that the inner scale is now 77, in 
agreement with our expectation. The exponentiation was 
achieved naturally in the present case. 

In assessing this result, we nee d to return to a deli- 
cate point in the derivation of Eq. ( 3.43| ) . The procedure 
involved summing all the terms of the order of unity, 
(powers of g ~ 1), w hile neglecting terms of O(e). How- 
ever, the sums (3.39-3.40) result in expressions containing 
(<? — 1) tx e. In other words, we end up with terms that 
appear of the same order as those neglec ted d uring the 
procedure. In order to justify the results (3.43) we must 
go back and analyze contributions of 0(e). These appear 
for example in contributions in which the "rungs" ap- 
pear on adjacent ladders, like case (b) in region II. The 
two "simple" rungs appearing in two adjacent ladders 
can be now considered as a single compounded rung. We 
focus now on the infinite set of diagrams in which this 
compounded rung repeats many times. The sum of such 
diagrams will again generate results containing (g — 1), 
and in the end will be responsible for terms of 0(e 2 ) in 
the scaling exponent. 

The reason for this phenomenon is the structure of the 
iterative solution. Sums of terms of order g have cancel- 
lations, leading eventually to a result of 0(e). The sum 
of terms of O(e) have a very similar structure, just we 
a redefined "rung" . Therefore they automatically gen- 
erate another factor of e by re-summation. This phe- 
nomenon repeats in higher orders, again by redefining 
what do mean by a "rung" . 

Thus Eq. ( |3~43| ) can be considered as the final result 
for the scaling of the fused correlation function of gradi- 
ent fields with the exponent correct to 0(e). In the next 
section we turn to the calculation of the scaling exponent 
of the unfused correlation function of the scalar field it- 
self. We will show that the exponents computed in both 
methods agree when the same objects are evaluated. This 
agreement is connected in the final section with the exis- 
tence fusion rules that control the asymptotic properties 
of unfused correlation functions when some coordinates 
are fused together. 



IV. ZERO MODES IN THE ANISOTROPIC 
SECTORS 

A. Calculation of the correlation functions 

In thi s section we consider the zero-modes of 
Eq. ( [2. 16 ). In other words we seek solutions Z n ({r m }) 



which in the inertial interval solve the homogeneous equa- 
tion 

n 

- iv)Vf Vf Z n ({r m }) = . (4.1) 

We allow anisotropy on the large scales. Since all the op- 
erators here are isotropic and the equation is linear, the 
solution space foliate into sectors {£, a} corresponding 
the the irreducible representations of the SO(d) symme- 
try group. Accordingly we write the wanted solution in 
the form 



Z n ({r m }) = J2 z n,eA{r m }) 



(4.2) 



where Z n ^ a is composed of functions which transform 
according to the (£, a)- irreducible representations of 
SO(d). Each of these components is now expanded in 
e. In other words, we write, in the notation of Ref. 



E„ 



eG nAr7 + 0{e 2 ) 



For e = Eq. (4.1) simplifies to 

n 



(4.3) 



(4.4) 



for a ny value of i, a. Next we expand the operator in 
Eq. (4.1) in e and collect the terms of 0(e): 



X]V?GrM>({r TO }) = V n E nAa ({r m }) , 



(4.5) 



where eV n is the first order term in the expansion of the 
operator in (4.1): 



V, 



a 
r jk r jk 



V?Vg , (4.6) 



where rjk = rj 

In solving Eq. (4.4) we are led by the following consid- 
erations: we want scale invariant solutions, which are 
powers of rjk- We want analytic solutions, and thus 
we are limited to polynomials. Finally we want solu- 
tions that involve all the n coordinates for the function 
E n 1 CT ; solutions with fewer coordinates do not contribute 



to the structure functions (2.E). To see this note that the 
structure function is a linear combination of correlation 
functions. This linear combination can be represented in 
terms of the difference operator Sj(r, r') defined by: 



5j{r,r')T({r m }) = F({r m })\ rj=r - T{{r rn })\ rj = T 



(4.7) 



Then, 
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S n (r 1 ,r 2 ...r r y n )=l[S ] (r J ,r' j )T({r m }) . (4.8) 



Accordingly, if T{{r m }) does not depend on r&, then 
Sk(Tk, r 'k)3~({r m }) — identically. Since the difference 
operators commute, we can have no contribution to the 
structure functions from parts of T that depend on less 
than n coordinates. Finally we want the minimal polyno- 
mial because higher order ones are negligible in the limit 
rjk <C A. Accordingly, E n ^^ with £ < n is a polynomial 
of order n. Consulting Appendix A for the irreducible 
representations of the SO(d) symmetry group, we can 
write the most general form of E n ^ tC! ^ up to an arbitrary 
factor, as 



(4.9) 



where [. . .] stands for all the terms that contain less than 
n coordinates; these do not appear in the structure func- 
tions, but maintain the translational invariance of our 
quantities. The appearance of the tensor B"Vo- Q '* °f Ap- 
pendix A is justified by the fact that E n ^^ must be sym- 
metric to permutations of any pair of coordinates on the 
one hand, and it has to belong to the £, a sector on the 
other hand. This requires the appearance of the fully 
symmetric tensor (A5). 

In light of Eqs. (4.5- 4~6|) we seek solution for 

Gn \{r m }) of the form 
G nM {{r m }) = ^Hi k tr ({r m })ln(r jk ) + H^({r m }) , 

j^k 

(4.10) 

where Hj k a ({r m }) and Hg ta ({r m }) are polynomials of de- 
gree n. The latter is fully symmetric in the coordinates. 
The former is symmetric in r^ , and separately in all 

the other {r m } m ^ hj 

Substituting Eq. ( 4.10| ) into Eq. ( |4.5| ) and collecting 
terms of the same type yields three equations: 



jk 
La 



V, ■ Vfci?2r. 



2[d-2 + r ifc -(V i -Vfc)]J3j' 



+ d-1 c 

i j^Lk 



2 K jk 



(4.11) 
(4.12) 



(4.13) 



Here Kj are polynomials of degree n — 2 which are sep- 
arately symmetric in the j, k coordinates and in all the 
other coordinates except j, k. In Ref. [|| it was proven 
that for £ = these equations possess a unique solution. 
The proof follows through unchanged for any £ ^ 0, and 
we thus proceed to finding the solution. 

By symmetry we can specialize the discussion to j = 1 , 
k = 2. In light of Eq. (4.12) we see that H\ 2 a must have 



at lea st a quadratic contribution in ri2- This guarantees 
that ( 4.10| ) is nonsingular in the limit r 12 — » 0. The only 
part of H\ 2 a that will contribute to structure functions 
must contain r$ . . . r n at least once. Since H} 2 a has to 
be a polynomial of degree n in the coordinates, it must 
be of the form 



rrl2 

n e,a — ' 12 ' 12 ' 3 



1,2 



(4.14) 



where [. . .]i j2 contains terms with higher powers of ri2 
and therefore do not contain some of the other coordi- 
nates rs . . .r n . Obviously such terms are unimportant for 
the structure functions. Since Hg 2 a has to be symmetric 
in ri, r 2 and r 3 . . .r n separately, and it has to belong to 
an £, a sector, we conclude that the constant tensor C 
has must have the same symmetry and to belong to the 
same sector. Consulting Appendix A, the most general 
form of C is 



° — at >2n.l.<j 



a l a 2 CK4 ■ ■ -02n 



2n-2,t,<r 
D 2n-4,,l,a 



(4.15) 



Substituting in Eq. (4.12) one finds 



(d + 2)Hj 



2d -2 



R ai...a 2 T, 



(4.16) 



I r «i„a2rttia2 7^1 
2'12'12° 



Substituting Eq. (4.14) and demanding that coefficients 



of the term 



will sum up to zero, we obtain 
2 



- 2(d+2)fl -2^r ' 

-2(d + 2)c=0; =^c = 



(4.17) 



The coefficient b is not determined from this equation due 
to possible contributions from the unknown last term. 
We determine the coefficient b from Eq. ( pi] ). After 
substituting the forms we find 



&2n 



\aB% 



■ ' 2n L"' J ^2rt,£,o- 

^271 DOl •••Q2n 

D 2n,l,a 



OL\Ct2 J^OL^a^ ...Q2 n 



J 1,2 



Recalling the identity (A6) we obtain 

, z n.i , 



Ad 



■[l-4o] 



2n-2,t,a J 

(4.18) 



(4.19) 



Finally we find that a is n, ^-independent, 

1 



2(d+2)(d- 1) 



(4.20) 



whereas b does depend on n and £, and we therefore de- 
note it as b n i 



bnj — 



(d + 1) 



4(d + 2)(d- 1) 



ZnJ ■ 



(4.21) 



In the next Subsection we compute from these results the 
scaling exponents in the sectors of the SO(<i) symmetry 
group with I < n. 
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B. The scaling exponents of the structure functions 

We now wish to show that the solution for the zero 
modes of the correlation functions T n (i.e Z n ) result in 
homogeneous structure functions S n . In every sector 



£ < n, a we compute the scaling exponents, and show 
that they are independent of a. Accordingly the scaling 
exponents are denoted Cn, and we compute them to first 
order in e. 



Using (4.7, 4.8), the structure function is given by: 



S^,«r(ri,ri; . . . ;r n ,r„) = A" 1 . --K^Kl 



(4.22) 



where A" ; = rf* — rf, and the function / is defined as: 



f^ir^r^Tj) ee [n -r^in-r^ ln|r 4 -r 3 \ + (n - rj)"* (n - fj)"* ln|r, ~r,| (4.23) 
-{n ~ (n - Tj)^ In \n - r 3 1 - {Ti - r 3 ) a * (r t - rj) a * In \r, -r 3 \ . 

The scaling exponent of S n j M can be found by multiplying all its coordinates by fi. A direct calculation yields: 

<5„,^ CT (/iri,/xri; . ..) = /J. n S n ,e, a {ri,ri; . . .) 

no i,j „o i,j 



2efx n ln M ]T A" 1 ■ ■ • AS- Af A?' [aB££°» + <^ B l-^\ + 0(0, 



^ A<r (ri,r i; . . .) - 2e^ n ln/iA? 1 . . . A«" ^[oB^ 01 - + Ki^K-iZ ] + °( e2 ) 



Using flAq), we find that 



£.<r > 



and therefore, we finally obtain: 

5 n (/iri,/iri; . . .) = - 2e[n(n - l)a + 6„^] hi/x}<S„(ri,Fi; . . .) + 0(e 2 ) 

= M C " > 5„(r 1 ,r 1 ;...) + 0(e 2 ) , 

The result of the scaling exponent is now evident: 



n(n — 1) 



(d+1) 



2(d + 2)(d-l) 4(d + 2)(d-l) 



z-aA + 0(e 2 ) 



n — e 



n(n + d) (d +!)£(£ + d - 2) 



2(d + 2) 2(d+2)(d-l) 



0(e 2 



For £ — this result coincides with [0]. This is the 
final result of this calculation. It is no tewort hy th at this 
result is in full agreement with (|3.43| ) and ( 3.21 ), even 
though the scaling exponents that appear in these result 
refer to different quantities. The way to understand this 
is the fusion rules that are discussed next. 



(4.24) 



C. Fusion Rules 



The fusion rules address the asymptotic properties of 
the fully unfused structure functions when two or more 
of the coordinates are approaching each other, whereas 
the rest of the coordinates remain separated by much 
larger scales. A full discussion of the fusion rules for the 
Navier-Stokes and the Kraichnan model can be found 
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in 0[LO|. In this section we wish to derive the fusion 
rules directly from the zero modes that were computed 
to O(e), in all the sectors of the symmetry group. In 
other words, we are after the dependence of the structure 
function <S n (ri, Fi; . . .) on its first p pairs of coordinates 
ri,ri; . . . ; r p , r p in the case where these points are very 
close to each other compared to their distance from the 
other n — p pairs of coordinates. Explicitly, we consider 
the case where ri , ri ; . . . ; r^, ^p-\-i , ^p+i ) • • • ] i" n , r n . 
(We have used here the property of translational invari- 
ance to put the center of mass of the first 2p coordinates 
at the origin). The calculation is presented in Appendix 
D, with the final result (to 0(e)) 

) (4-25) 

v 

= y ] ^ ] '0j,CT / 'S'p,j.(T'( r ii r ii---i r jn r p) ■ 



In this expression the quantity tpj a , is a function of all 
the coordinates that remain separated by large distances, 
and 



Jt, 



= max{0, p + £ — n} , t < n 



(4.26) 



We have shown that the LHS has a homogeneity expo- 
nent Cn ] - The RHS is a product of functions with homo- 
geneity exponents £p and the functions if) a a i. Using the 
linear independence of the functions Sp_j^> we conclude 
that tpj a , must have homogeneity exponent £„ — £p . 
This is precisely the prediction of the fusion rules, but 
in each sector separately. One should stress the intu- 
itive meaning of the fusion rules. The result shows that 
when p coordinates approach each other, the homogene- 
ity exponent corresponding to these coordinates becomes 
(i) 

simply Q as if we were considering a p-order correlation 
function. The meaning of this result is that p field am- 
plitudes measured at p close-by coordinates in the pres- 
ence of n — p field amplitudes determined far away be- 
have scaling-wise like p field amplitudes in the presence 
of anisotropic boundary conditions, closing we note that 
the tensor functions ijtj a i do not necessarily belong to 
the j, a' sector of SO(d). 



V. SUMMARY AND DISCUSSION 

One of our motivations in this paper was to under- 
stand the scaling properties of the statistical objects un- 
der anisotropic boundary conditions. The scaling expo- 
nents were found for all I < n, cf. Eq.(l.l). We found 
a discrete and strictly increasing spectrum of exponents 
as a function of i. This means that for higher £ the 
anisotropic contributions to the statistical objects decay 
faster upon decreasing scales. In other words, the statis- 
tical objects tend towards locally isotropic statistics upon 
decreasing the scale. The rate of isotropization is deter- 
mined by the difference between the £ dependent scaling 



exponents, and is of course a power law. The result shows 
that the ^-dependent part is n-independent. This means 
that the rate of isotropization of all the moments of the 
distribution function of field differences across a given 
scale is the same. This is a demonstration of the fact 
that the distributions function itself tends towards a lo- 
cally isotropic distribution function at the same rate. We 
note in passing that to first order in e the £ dependent 
part is also identical for £2, a quantity whose isotropic 
value is not anomalous. For all I > 1 also is anoma- 
lous, and in agreement with the n — 1 value of Eq.(l.l). 
Significantly, for £2 we have a nonperturbative result that 
was derived in j?]], namely 



(i) 



1 r 
2 



(5.1) 



{2-d 



4(d + e- \)£{d + £-2) 



1 



valid for all values of e in the interval (0,2) and for all 
£ > 2. Th is exact result agrees after expanding to 0(e) 
with (|3i|) for n = 1 and £ = 2 . 

Our second motivation was to expose the correspon- 
dence between the scaling exponents of the zero modes 
in the inertial interval and the corresponding scaling ex- 
ponents of the gradient fields. The latter do not depend 
on any inertial scales, and the exponent appears in the 

combination (A/£) ™ where £ is the appropriate ultravi- 
olet inner cutoff, either A or 77, depending on the limiting 
process. We found exact agreement with the exponents of 
the zero modes in all the sectors of the symmetry group 
and for all values of n. The deep reason behind this 
agreement is the linearity of the fundamental equation of 
the passive scalar (2.1). This translates to the fact that 
the viscous cutoff rj, Eq.(2.12) is n and £ independent, 



and also does not depend on the inertial separations in 
the unfused correlation functions. This point has been 
discussed in detail in |f^JTl[] . In the case of Navier-Stokes 
statistics we expect this "trivial" correspondence to fail, 
but nevertheless the "bridge relations" that connect these 
two families of exponents has been presented in Jl0[ for 
the isotropic sector. Finally we note that in the present 
case we have displayed the fusion rules in all the £ sec- 
tors, using the 0(e) explicit form of the zero modes. We 
expect the fusion rules to have a nonperturbative valid- 
ity for any value of e. It would be interesting to explore 
similar results for the Navier-Stokes case. 
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APPENDIX A: ANISOTROPY IN 
D-DIMENSIONS 

To deal with anisotropy in d-dimensions we need clas- 
sify the irreducible representations of the group of all 
(i-dimensional rotations, SO(d) jl^], and then to find a 
proper basis for these representations. The main linear 
space that we work in (the carrier space) is the space of 
constant tensors with n indices. This space possesses a 
natural representation of SO(d), given by the well known 
transformation of tensors under d-dimcnsional rotation. 

The traditional method to find a basis for the irre- 
ducible representations of SO(<i) in this space, is using 
the Young tableaux machinery on the subspace of trace- 
less tensors jl2| . It turns out that in the context of the 
present paper, we do not need the explicit structure of 
these tensors. Instead, all that matters are some relations 
among them. A convenient way to derive these relations 
is to construct the basis tensors from functions on the 
unit d-dimensional sphere which belong to a specific ir- 
reducible representation. Here also, the explicit form of 
these functions in unimportant. All that matters for our 
calculations is the action of the Laplacian operator on 
these functions. 

Let us therefore consider first the space Sd of functions 
over the unit d-dimensional sphere. The representation 
of SO(d) over this space is naturally defined by: 



(Al) 



where 9(u) is any function on the d-dimensional sphere, 
and 1Z is a rf-dimensional rotation. 

Sd can b e sp anned by polynomials of the unit vector u. 
Obviously (Al) does not change the degree of a polyno- 
mial, and therefore each irreducible representation in this 
space can be characterized by an integer £ = 0,1,2,..., 
specifying the degree of the polynomials that span this 
representation. At this point, we cannot rule out the 
possibility that some other integers are needed to fully 
specify all irreducible representations in Sd and therefore 
we will need below another set of indices to complete the 
specification. 

We can now choose a basis of polynomials {Y^ jCT (u)} 
that span all the irreducible representations of SO(d) over 
Sd- The index a counts all integers other than £ needed 
to fully specify all irreducible representations, and in ad- 
dition, it labels the different functions within each irre- 
ducible representation. 

Let us demonstrate this construction in two and three 
dimensions. In two dimensions a is unneeded since all 
the irreducible representation are one-dimensional and 
are spanned by Yi(u) = e u< ^ with <f> being the angle 
between u and the the vector e\ = (1,0). Any rota- 
tion of the coordinates in an angle <po results in a mul- 



tiplicative factor e 4 ^ . It is clear that Yi(u) is a poly- 
nomial in u since Yi{u) = [u ■ p] where p = (l,i). In 
three dimensions a — m where m takes on 2£ + 1 values 
m = -£,-e+l,...,e. Here Y^ m cx e im ^Pp(cos 9) where 
4> and 9 are the usual spherical coordinates, and P™ is 
the associated Legendre polynomial of degree £ — m. Ob- 
viously we again have a polynomial in u of degree £. 

We now wish to calculate the action of the Laplacian 
operator with respect to u on the Yi a (u). We prove the 
following identity: 



d a d a Y^ a {u) = -£{£ + d - 2)Y e ^(u) 



(A2) 



One can easily check that for d — 3 (|A2|) gives the fac- 
tor £{£ + 1), well known from the theory of angular- 
momentum in Quantum Mechanics. To prove this iden- 
tity for any d, note that 



e d 2 \u\% a (u))=Q 



(A3) 



This follows from the fact that the Laplacian is an 
isotropic operator, and therefore is diagonal in the Yi^. 



The same is true for the operator \u\ d 



But this 



operator results in a polynomial in u of degree £ — 2, 
which is spanned by Ygj . CT / such that £' < £ — 2. Therefore 
the RHS of (|A3|) must vanish. Accordingly we write 



0VPM<i) + 2d c V|cr^, (T + \u l \d'Y^(u) = o . 

(A4) 

The second term vanishes since it contains a radial 
derivative u a d a operating on Y^ )Cr (u) which depends on 
u only. The first an d th ird terms, upon elementary ma- 
nipulations, lead to ( A2). 

Having the Y^^(u) we can now construct the irre- 
ducible representations in the space of constant tensors. 
The method is based on acting on the Yg j(7 (u) with the 
isotropic operators u a , d a and 8 al3 . Due to the isotropy 
of the above operators, the behavior of the resulting ex- 
pressions under rotations is similar to the behavior of the 
scalar function we started with. For example, the tensor 
fields 5 a ^Yi il7 (u), d a d^Yi^(u) transform under rotations 
according to the (£, cr) sector of SO(d). 

Next, we wish to find the basis for the irreducible rep- 
resentations of the space of constant and fully symmetric 
tensors with n indices. We form the basis 



D"l,-, 

D t,o,n 



d ai ...d an u n Yt !<T (u), £<n 



(A5) 



Note that when £ and n are even (as is the case invari- 
ably in this paper), B^'™'°" no longer depends on u, 
and is indeed fully symmetric by construction. Simple 
arguments can also prove that this basis is indeed com- 
plete, and spans all fully symmetric tensors with n in- 
dices. Other examples of this procedure for the other 
spaces are presented directly in the text. 

Finally let us introduce two identities involving the 
B n ,e,a which are used over and over through the paper. 
The first one is 
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= Zn*B%~% , (A6) 
z„. £ = [n(n + d-2)-i{l + d-2)\ . (A7) 



It is straightforward to derive this identity using ([Aj). 
The second identity is 



{a m },mjti,j 
e,a,n-2 



I < n - 2 . (A8) 



This identity is proven by writing u n in ( A5) as 



and operating with the derivative on u 2 . The term ob- 
tained as u 2 d ai . . .d an u n ~ 2 Yi^(u) vanishes because we 
have n derivatives on a polynomial of degree n — 2. It 
is worthwhile noticing that these identities connect ten- 
sors from two different spaces. The space of tensors with 
11 indices and the space of tensors with n — 2 indices. 
Nevertheless, in both spaces, the tensors belong to the 
same (I, a) sector of the SO(d) group. This is due to the 
isotropy of the contraction with S aia2 in the first identity, 
and the contraction with S ai0lj in the second identity. 



APPENDIX B: PROOF OF EQ. (3.32) 



contributes a term proportional to gAffl if the p vector 
associated with this rung is the largest among all the 
p vectors associated with rungs appearing to the right 
of it. Otherwise the contribution is proportional to g. 
The weight of the contribution is obtained as a factor 
c < 1 which reflects the proportional fraction of the vol- 
ume of (p 1 ,p 2 ■ ■ .)-space in which the associated ordering 
of the p vectors is valid. For example, if the rung with 
the largest p vector is in the extreme right, then all the 
other rungs contribute terms proportional to g. Thus a 
diagram with s rungs ordered in this manner contributes 
with a weight C = 1/s. Therefore 



Di, a = - 
s 



(CI) 



The recurrence relation for D m , s with m > 1 is derived 
by inserting an additional rung which is associated with 
the largest p vector in any one of the (s + 1) possible po- 
sitions available in a diagram with s rungs. After some 
combinatorial calculations of the weights one finds 



1 s_1 



q,m • 



(C2) 



In case a of region II where Pi > P2 > k s the analytic 
expression for A 2 can be simplified to 



A' 



1 



x 2.0i0i - 2 

Using the identities 



Pi Pi 



d V2 n/3j/3 . 



dp = Q(d) , n d = Q(d)- 



dpp a f = s a(3 n(d) , 



d-i 



we compute 



(Bl) 

(B2) 
(B3) 

(B4) 



Substituting this form into the double rung ladder dia- 
gram results, after contracting all the indices of Ai and 
A 2 , in a form identical to Eq. (3.16) for A in the one rung 
ladd er dia gram. This leads directly to the final equation 
Eq. (|3~32l) . 



APPENDIX C: DOUBLE RESUMMATION 
1. Calculation of D m ,s 

In this Appendix w e disc uss the calculation of the co- 
efficients D m ^ s in Eq. (3.35), and the actual resummation 
of that equation. 

Firstly we need to introduce rules to evaluate the rungs 
in the general ladder diagram that appears in the ex- 
pansion. The rule is actually quite simple: every rung 



Together with (CI) this gives 

s-l 



D 



-E 1 

S ^2 I 2 £i * 

1 x - 1 q ^-^ 1 v - ^ J 

s ^ q3 ^ <72 , qi 

93=3 92=2 9i = l 

The general structure of D m s now becomes obvious 



D 



s " q_2 " qi 

92=2 1 91 = 1 
- s-l ^ 93-1 g 2 -l j 



etc. 



(C3) 
(C4) 
(C5) 



2. Higher order terms in A 



Consider the equation (|3.38| ) 

00 ~ e s — 1 



K2(9,A)=A^ilj:l 



Observing that 

1 1 



<q 2V^\s-q q 



S =2 3 9=1 q 



1 1 



(C6) 



and 

s \q s-q 
we end up with 



q(s - q) 



(C7) 



(C8) 
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CO oo 



91=1 g 2 = l 



91+92 



9i 92 



(C9) 



where we relabeled q — > q\ and s — q —> q2 and changed 
correspondingly the limit of summation over q2- Thus 



K 2 (~g,A) = ^ 2 [E-] 2 = \[-AHl-9) 

9=1 ^ 

The terms proportional to A 3 give 

■o -„ s-1 1 92-1 



71 " g 2 " 91 

s=3 92=2 9i = l 



(CIO) 



(Cll) 



We can rearrange the sums by summing over q% — 
n — q\ — 172 instead of n. Using relationships similar to 
flC7|) and Q we find 



v X1 ^91+92+93 



-4 va; v-^j 

^ = ^ 2 EEE 

91=1 92=1 93=1 

Obviously this leads to 

9 



919293 



(C12) 



^^) = ^[Ea 3 .i f [-Aln(l- 

9=1 



(C13) 



The general structure is now clear, leading to Eq. (3.39) 



APPENDIX D: DERIVATION OF THE FUSION 
RULES 



In this appendix we derive the fusion rules ( 4.25 ) . Con- 
sider a fully unfused structure function with n coordi- 
nates, such that p of them are separated from each other 
by a typical distance r, whereas n — p coordinates are 
separated from them and from each other by a typical 
distance R, and R^$> r. We want to compute the asymp- 
totic properties o f S n and show that to leading order in 
r/R we find Eq.(4.25). For homogeneous ensembles we 
can shift the origin to the center of mass of the p coor- 
dinates. In this case we have Tj <C ri for every j < p 
and i > p. Our aim is to separate the dependence on the 
small distances from the dependence on the large dis- 
tances. We will see that some of the terms in S n lend 
themselves naturally to such a sep arati on, and some call 
for more work. We start from Eq.(4.22), and compute to 
first order in r/R: 



r^ir^r^fj) ee ( n - rj r(n - rj ) a ' In |r, - Tj \ + (n - r^in - ln|r 4 - r,| 
-(n - rj) a ' {n - rj) a * In |r; -fj\ - (r l - rj) " (r* - r,-)^ In |Ej - r, 



^"^(r^r^A 



(4.1) 



and so, if rj, rj <§C r^, for j = 1, . . . ,p , i = p + 1, . . . , n then the first order in e of S n will contain 3 types of terms: 



no t,j 



no t,j 



E . . . A«» r^(r, l r <1 r il r j )K^ + 6„,M Q ^ 

i<»#j'<p 

no i,j no ij 



E Ai 1 • • ■ Ap p /<"°« (r is r,, r^HaS^ + ] A^ 1 ■ • • A- 



no i,j no i,j 

E Af r a n^v u v 3 ,r 3 )[aB^ +b nA 5 a * a * B^~g] 

l<j<p,p<i<n 

no i no 



A a i A ft p A a p+i a^^ a,/3 / rA\aB ai ;"P"' an 4- ft „f^>i- Q n 



l^J ' <p,p<z<n 



no 2 j 



J 3 = V Af...A^/^(r i! r i ,r 3 .,f j )[aB, 



p<2 j<n 



no 

+ b nlO tin-It., 
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A Qp 



E 

p<i,j<n 



^p+l 



We note that of these three terms only 1% has a nontrivial mixing of small and large coordinates, and indeed it is the 
only term in which the expansion (4.1) was employed. Collecting terms we find 



where: 



no z j 



6ta(ri,ri; . . . ; r„,f„) = A" 1 . . . A^B^- a - 

no i,j 

+e^2 r ai (r< 1 fi,rj,r J )[a^ 1 '" a » + 6 B ,ii° ,a * +eA? 1 . . . A^q*™ * + 0(e 2 ), 



5oi...a, _ a^p+i a a n nai...a„ _ \^ „ R ai...a p 

#p p - A P +i •■•V i, n ,v _ Z^'-'-- T 'V.."' 

j-max{0,]3+£-n} <7 ; 
P-2 



D a 3-«p A a P+l A Q " R«3-«I1 _ 

-°p-2 ~~ ^p+l ■ ' ■ D 



j-raax{0,p+£-n} cr' 



P 



£ A^ 1 .. . A£» r^r.f.r^JK^" + & n ,,<^ 

p<ij'<n 

no i no ij 

+ P+l ••• ZA n 5 " 3 I 1 *' I i)[ ar> n,e,<r + V - D n-2/, CT J 

l<j<p,p<i<n 

E E-,-'w • 



In these expressions we use the fact B p is a fully symmetric tensor and therefore can be again expanded in terms of 
the basis functions B pjj(T / with coefficients that depend on the large separations. The sums on the right hand sides 
run between j = max{0, p + t — n} and p because not all the basis functions can appear when p + 1 — n > 0. This can 
be checked by contracting the B p with (n — £ + 2)/2 delta function. This contraction vanishes since it contains a factor 
Zij.. On the other hand the contraction results in a tensor with p + 1 — n — 2 indices, and therefore all corresponding 
coefficients of j < p + £ — n — 2 must vanish. To proceed we establish the following identity: 

M c . ff , =d . )CT/ (4.2) 
Zn,t 

The identity is proven by the following calculations: 

a;^ 1 . . . K n K:c n = d ai --- d^vp u -n^ ■ ■ ■ K n )9 ap+1 . . . ^-r^^ 

j.a' 



A;^ 1 . . . K n K-2Z = d ° a ■ ■ ■ d^u^u-^^lX 1 . . . A^}d ap+1 . . . d Qn u n - 2 Y La (u) 

= as ... 8 a "U P dj.a'Yj.v, (u) . 

Denote now 

f(u) = u-'PiA^ 1 . . . A^]d p+1 . . . d n u n Y^(u) = Cj,<r>Yi,«' («) 
g(u) = u-^lA^X 1 . . . A^]d p+1 . . . d n u n - 2 Ye,„(u) = £ d^Y^(u) 



1G 



To obtain (4.2), operate with u 2 d 2 on f(u). On one hand, we get: 

,,2 



d 2 f(u) = ~3(J + d - 2 )^Y jt a'(u) 



but on the other hand, we have: 



u 2 d 2 f(u) = u 2 d*[u- p [A^ + + 1 . . . A^]d p+1 . . . d n u n Y it(T (i 
= -p(-p + d- 2)f(u) - 2p 2 f{u) + z n/ g(u) 
= -pip + d- 2)f(u) + z ni eg(u) . 

Equating the two expressions, and projecting over the (J, a') sector, we obtain: 

+ d- 2)c 7> < = -p{p + d- 2)c jt<7 > + z n ,idj,v> 

\p(p + d-2)-j(j + 2-2)} z P: 



Zn,l 



Recalling Eq.(4.21), b p ,j = b ni iz Pj j / z n> i and we may write to leading order in r/R: 



S n (ri,f i; ...;r„,r„) = £ ^\ A',' . . . A" KX^ 

no i,j no i,j 

+£ ^A^r ) (r I ,r ! ,r J ,r 3 )[ fl B^+5 PJ ^ a] 
+ E £ ■ ■ ■ + 0(e 2 ) 

j-max{0,p+f-Ti} u' 
P 

E J2(c 3 ,*> +£e^ a ,)Sp, 3 , <T >(ri,f 1 ;...;r p ,r p ) + 0{e 2 ) 

j-max{O t p+£-n} tj' 



From this follows Eq.(4.25) 
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